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In the paper, the thermodynamic properties of the CaLi2 superconductor have been studied. 
It has been assumed the high value of the pressure: p = 100 GPa. The thermodynamic critical 
field (He) and the specific heats for the superconducting and normal state (C s and C N ) have 
been obtained. The calculations have been made within the framework of the Eliashberg approach 
for the wide range of the Coulomb pseudopotential: fi* G (0.1,0.3). It has been shown that the 
low-temperature critical field and the specific heat jump at Tc strongly decrease if fi* increases. 
The value of the dimensionless ratio Rc = (C s — C N ) /C N at the critical temperature differs 
significantly from the value predicted by the BCS model. In particular, the parameter Rc decreases 
from 1.83 to 1.73 with the increase of the Coulomb pseudopotential value. In the paper, the 
interpolation formulas for the functions Tc ((J,*) and Rc (n*) have been also given. 
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Calcium under the compression exhibits the anomalous 
behavior associated with the structural transitions. In 
particular, Ca transforms from the fee to bec structure 
at 19.5 GPa [l[. The simple cubic structure exist in the 
pressure (p) range from 32 GPa to 113 GPa. Above exist 
the Ca-IV phase Q ■ Several other higher pressure phases 
like Ca-V (at p = 139 GPa @), Ca-VI (at p = 158 GPa 
0) and Ca-VII (at p = 210 GPa [|) have been also 
confirmed into the experimental way. 

At p = 44 GPa (the sc phase), the superconductiv- 
ity with the critical temperature (Tc) close to 2 K has 
been observed The critical temperature in calcium 
increases with the pressure and reaches the maximum 
value equal to 25 K at 161 GPa [(| - the highest value of 
the critical temperature for simple elemental supercon- 
ductors. 

Calcium can create the compounds with many other 
elements. The theoretical studies have predicted that the 
hydrogen-rich compound like CaHg (p — 150 GPa) can be 
the high temperature superconductor with the extremely 
high value of the critical temperature: Tc <G (220 — 235) 
K However, this result is still not confirmed by the 
experimental data. 

On the other hand, the superconducting state in CaLi2 
has been experimentally observed at the pressure above 
11 GPa. From 11 GPa to 40 GPa, the critical tempera- 
ture increases to maximum value equal to 13 K, then Tc 
slowly decreases to 9 K at 81 GPa 

The theoretical works suggest that CaLi2 compound 
above 81 GPa can also possesses the interesting super- 



conducting properties. In particular, the critical temper- 
ature assumes the high value @. Therefore in the pre- 
sented paper, we have study the thermodynamic proper- 
ties of the superconducting state of CaLi2 at the pressure 
of 100 GPa. 

In the considered case, the electron-phonon coupling 
constant is large: A = 2.73 @; so the mean-field BCS 
theory can not be used [l(J. For this reason, the calcu- 
lation have been conducted within the framework of the 
Eliashberg formalism [ll|. In particular, the thermody- 
namic critical field and the specific heats for the normal 
and superconducting state have been determined. 

The Eliashberg equations on the imaginary axis can be 
written in the following form: 
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where the symbols A„ = A (ioj n ) and Z n = Z (iui n ) de- 
note the order parameter and the wave function renor- 
malization factor, respectively. The n-th Matsubara fre- 
quency is given by the expression: uj n = Jj (2n — 1), 
where j3 = 1/fcsT, and fc^ is the Boltzmann constant. 

The electron-phonon pairing kernel K (n, m) is given 
by: 
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where the Eliashberg function (a 2 F (O)) and the maxi- 
mum phonon frequency (r2 max = 113.1 meV) has been 
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FIG. 1: The order parameter and wave function renormalization factor on the imaginary axis for the selected values of the 
temperature and the Coulomb pseudopotential. In the figures (A)-(C), the first 150 values of the order parameter have been 
presented; (D)-(F) the 150 values of the wave function renormalization factor. 



determined in the paper In particular, the crystal 
structure for CaLi2 compound has been obtained by us- 
ing the ab initio evolutionary algorithm USPEX [l2| , [l3| , 
14 . This method has been also in detail discussed in the 
papers [l5[ , [l6| , [13] ■ The structural relaxations and the 
electronic properties calculations have been carried out 
using the density functional theory [IH, [lj|, [2(J (the 
VASP code [Iljj). The lattice dynamics and the electron- 
phonon coupling have been analyzed with help of the 
QUANTUM-ESPRESSO package [H. 

In the presented paper, the depairing electronic inter- 
action has been described by the Coulomb pseudopoten- 
tial (//*); the symbol 9 denotes Heaviside unit function 
and uj c is the cut-off frequency (w c = 3f2 max ). 

Due to the absence of the experimental value of the 
critical temperature, we assume the wide range of the 
Coulomb pseudopotential: fi* £ (0.1,0.3). This choice is 
connected with the fact that the value of the Coulomb 
pseudopotential in Ca and Li is high [23|, (24|. For ex- 
ample: [/iciSlGPa = °' 24 and [Mc4=22.3GPa = °- 22 OT 

[Mc4=k 7G p a = 0-36 MM- 

The Eliashbcrg equations have been solved for 2201 
Matsubara frequencies b y u sing of the iteration method, 
described in the papers [27| and [28j |. In the considered 
case, the functions A„ and Z n are stable for T > Tq, 
where the temperature To is equal to 2.32 K. 

In Fig.[TJ the solutions of the Eliashberg equations on 
the imaginary axis have been presented. We have chosen 
the temperature range from To to Tq for the selected 
values of Coulomb pseudopotential. 



In Fig.Q] (A-C), the dependence of the order parame- 
ter's values on the successive Matsubara frequencies has 
been plotted. It is easy to see that the maximum value 
of the order parameter (A m= i) is decreasing with the 
growth of the temperature; the half-width of the A m 
function becomes successively smaller. The last prop- 
erty suggests that together with the growth of the tem- 
perature, less of the Matsubara frequencies give relevant 
contribution to the Eliashberg equations. 

In Fig.[T](D-F), the wave function renormalization fac- 
tor on the imaginary axis has been shown. In contrast 
to the maximum value of the order parameter, Z m= \ 
slightly increases if the temperature increases. The half- 
width of the Z m function decreases with T . However, 
in comparison with the order parameter, this decrease is 
much slower. 

We notice that the increase of the Coulomb pseudopo- 
tential strongly lowers the values of the order parameter 
and slightly influences on the wave function renormaliza- 
tion factor. 

The thermodynamic properties of the superconducting 
phase depend on the free energy difference between the 
superconducting and the normal state (AF = F s — F N ) 
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where p (0) represents the value of the electron density 
of states at the Fermi surface. The symbols Z% and Z^ 
denote the wave function renormalization factors for the 
superconducting and the normal state, respectively. 

In Fig.[3J the dependence of AF/p(0) on the temper- 
ature for selected values of the Coulomb pseudopotential 
has been plotted. The negative values of the considered 
ratio prove that the superconducting state is thcrmody- 
namically stable to the critical temperature. Addition- 
ally, it can be noticed that AF (To) strongly decreases 
with (ji*: 



[AF (To)U =0 .3 / [AF (T )]^ =0A ~ 0.41. 



(5) 
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FIG. 2: The free energy difference between the superconduct- 
ing and the normal state as a function of the temperature for 
selected values of the Coulomb pseudopotential. 

The thermodynamic critical field can be calculated in 
the following way [3(| : 



V-8tt [AF/p l 



(6) 



In Figjjp the temperature dependence of the ratio 
He I v P (0) has been shown. It is easy to see that the 
low temperature value of the critical field strongly de- 
creases if the Coulomb pseudopotential increases: 



[Hc(T )]^ , 3 /[H c (T )]^ : 



0.65. 



(7) 



The specific heat difference between the superconduct- 
ing and the normal state (AC = C s — C N ) should be 
calculated by using the expression: 



AC (T) _ 1 d 2 [AF/p (0)] 
W(0) ~ ~P d(k B T) 2 



(8) 
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FIG. 3: The thermodynamic critical field as a function of the 
temperature for selected values of the Coulomb pseudopoten- 
tial. 



Whereas, the specific heat in the normal state can be 
obtained from the formula: -, — W- = where the Som- 
merfcld constant is given by: 7 = |7r 2 (1 + A). 



C s / k B p(0) meV 

C N /k B p(0)meV o0 o°°° ; 



- 1 H*=0.2 
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FIG. 4: The dependences of the specific heats in the super- 
conducting and the normal states on the temperature for se- 
lected values of the Coulomb pseudopotential. The vertical 
line indicates the position of specific heat's jump at Tc. 

In Fig. |4l the form of the specific heats has been shown. 
Let us notice that at low temperatures, the specific heat 
in the superconducting state is exponentially suppressed. 
For the higher temperature, C s rapidly increases and 
much exceed the value of C N . At the critical temperature 
the characteristic jump has been observed. 

Furthermore, we can see that at the critical tempera- 
ture the specific heat's jump decreases with the growth 
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of the Coulomb pscudopotcntial: 



[&C(Tc)l* =Q . 3 /[AC(T c )}^ 0A 



0.59. 



(9) 



The knowledge of the functions C s and C N enables 
the determination of the parameter: 



Rc 



AC (7b) 
C N (T c ) ■ 



(10) 



For CaLi 2 compound under the pressure at 100 GPa, 
we have obtained: R c G (1.83,1.73), if ^* G (0.1,0.3). 
The above result strongly differs from the value predicted 
by the BCS model: [Rc] BCS = 1.43 0. The observed 
difference between the prediction of the BCS and Eliash- 
berg theory is connected with the existence of the strong- 
coupling and retardation effects in the CaLi2 compound. 
These effects are measured by the parameter: /csTc/win, 
where the symbol aj\ n denotes the logarithmic phonon fre- 
quency [13 ■ In the BCS limit, we have: ksTc/^in — > 0. 
On the other hand, in the framework of the Eliashbcrg 
theory, we can obtain: ksTc/^Xix G (0.656,0.404) for 
/j* G (0.1,0.3). We notice that the critical temperature 
has been calculated from the formula determined in the 
Appendix and u!\ n = 3.87 mcV. 

In order to determine the value of Rc, it has been 
necessary to made the complicated calculations. For this 
reason, we present the convenient interpolation formula 
also for Rc in order to reconstruct the precise numerical 
result: 
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In summary, we have used the Eliashbcrg equations to 
calculate the thermodynamic critical held and the spe- 
cific heats of CaLi2 superconductor under the pressure 
at 100 GPa. We have shown that the low-temperature 
value of the critical field and the value of the specific 
heat jump at Tc strongly deceases if the Coulomb pseu- 
dopotential increases. The ratio Rc for the characteristic 
values of the functions C s and C N can not be described 
correctly by the conventional BCS theory In particular, 
R c G (1.83,1.73), where fx* G (0.1,0.3). The numerical 
results for the functions Tc (^*) and Rc (fJ*) have been 
supplemented by the interpolation formulas. 

Acknowledgments 

The calculations have been based on the Eliashberg 
function sent to us by: Prof. Yanming Ma and Ph.D. Yu 
Xie to whom we are very thankful. 

The authors are grateful to the Czestochowa University 
of Technology - MSK CzcstMAN for granting access to 
the computing infrastructure built in the project No. 
POIG.02.03.00-00-028/08 "PLATON - Science Services 
Platform" . 

Appendix 

The interpolation formula for the critical temperature 
takes the form: 



k R T, 



B±C 



/l/2 L45 CXP 



-1.172(1 + A) 



(Al) 



where /i and f% denote the strong-coupling correction 
function and the shape correction function, respectively 
HI: 



aT 



and, / 2 



1 A 2 



The quantities Ai and A 2 have the form: 



A 2 + A2 
(Al) 



FIG. 5: The value of the ratio Rc as a functions of the 
Coulomb pseudopotential. The circles represent the exact 
Eliashberg result. The solid line has been obtained on the 
basis of the expression (|11[) . The BCS result has been also 
shown. 

In Fig [H we have presented the ratio Rc as a function 
of the Coulomb pseudopotential. It is easy to see that the 
formula (jTTJ) reproduces correctly the numerical data. 



A 1= 2-2.6/x* and, A 2 = 0.31 + 1.4 1// [*—=- . 

V w ln / 

(Al) 

The parameter ,/uJ2 is equal to 34.95 meV. 

The formula (|A1|) has been prepared basing on 300 
exact values of the function Tc (/a*). We have used the 
last square method and we have considered the range of 
fi* from 0.1 to 0.3. 
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